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We consider modified f(R) gravity which may unify R m early-time inflation with late-time ACDM 
epoch. It is shown that such model passes the local tests (Newton law, stability of Earth-like 
gravitational solution, very heavy mass for additional scalar degree of freedom) and suggests the 

£ — | realistic alternative for General Relativity. Various scenarios for future evolution of f(R) ACDM 

f^*) , era are discussed. 
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I. INTRODUCTION 



Modified gravity is considered as a very interesting alternative proposal for dark energy. The attractive property 
of gravitational dark energy (for a review, see [!}) is the fact that one should not introduce some strange matter 
with negative pressure to describe the late-time cosmic acceleration. The (sudden) change of decelerated expansion 
to accelereted one is explained by the change of the properties of a gravitational theory in the course of the universe 
evolution. In other words, some sub-leading gravitational action terms [H, 0, E| m ay become essential ones at the 
tIh ■ late universe. That is why there is much activity in the study of different versions of modified f(R) gravity with 
applications to dark energy cosmology [H, H, H, 0) IE S 0> H] • 
! Recently, a very realistic modified f(R) gravity which evades the solar system tests was proposed in ref.@ (for 
J> 1 related discussion, see [10]). On the same time, such theory leads to an effective ACDM epoch which complies 
with observational data with the same accuracy as the usual General Relativity (GR) with cosmological constant. 
Generalization of the model @ done in ref. [ll[ suggests a quite natural unified description of early-time inflation with 
late-time acceleration following the earlier proposal of ref. 3 . In the present work we propose another class of modified 
f(R) gravity which unifies R m inflation with ACDM era. It is shown that such theory may pass the local tests. The 
cosmological properties of such model are studied. Some speculative remarks about the possibility of a (future/past) 
anti-gravity phase are made. It is also shown that such theory may be more friendly with future observational data: 
unlike to usual GR the current ACDM epoch induced by such theory may enter to future quintessence/ transient 
phantom era by the effective reconstruction of the action. Alternatively, one can have the eternal ACDM epoch as a 
!>" ■ result of such reconstruction. 

X: 

C<3 1 II- UNIFYING R m INFLATION WITH ACDM COSMOLOGY 

We start from the following action of general f(R) gravity 

S=± J ' d A x^—g(R + f(R)) . (1) 
The equation of motion in /(i?)-gravity with matter is given by 

\g^F{R) - R^F'(R) - g^UF'(R) + V^J U F'{R) = -yT (m)M „ . (2) 
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Here F(R) = R + f(R) and TV m )^„ is the matter energy-momentum tensor. By introducing the auxilliary field A one 
may rewrite the action ([1]) in the following form: 

S = J_ / d 4 xV^{(l + f'(A)) (R-A)+A + f(A)} . (3) 



As is clear from ([3]), if F'(R) = l + f'(R) < 0, n 2 s = n 2 /F'(A) becomes negative and the theory enters the anti-gravity 
regime. Note that it is not the case for usual GR. 

Recently a viable £(R) model has been proposed in ref.[9((for other recent proposals/study of properties of viable 
f(R) gravity, see @,B In this model f( R ) 

is chosen to be 

m 2 ci (R/m 2 ) n 

JHS{R) = , p/ 2 yn . , , (4) 

c 2 (i?/ m ) +1 



which satisfies the conditions 



lim Jhs(R) = const , lim fns(R) = . (5) 

it— >oo it^O 



The second condition means that there is a flat spacetime solution (vanishing cosmological constant). The estimation 
of ref.Q suggests that R/m 2 is not so small but rather large even in the present universe and R/m 2 ~ 41. Hence, 

f / m m2c i , m2c i ( R Y n /-n 
fHs(R) — + —3-32 . ( 6 ) 



C2 

which gives an "effective" cosmological constant — m 2 ci/c 2 and generates the late-time accelerating expansion. One 
can show that 

ff2 ^ m 2 ^ ^ (70km/g . pc) 2 _ ( 10 -33 e y)2 _ (?) 
C2 

Hence, the above model describes an effective ACDM cosmology. 

Although the model [9( is very succesful, early time inflation is not included there. We have suggested the modified 
gravity model to treat the inflation and the late-time accelerating expansion in a unified way [ll| . We have considered 
a simple extention of the model @ to include the inflation at the early universe. In order to generate the inflation, 
one may require 

lim f(R) = -A, . (8) 

R — >oo 

Here A^ is an effective cosmological constant at the early universe and therefore it is natural to assume Aj ^> 
(lO~ 33 eV) . For instance, it could be Aj ~ xo 20oj38 (eV) 2 . In order that the current cosmic acceleration could be 
generated, let us consider that currently f(R) is a small constant, that is, 

f(R ) = -2i? , f'(Ro) ~ . (9) 

Here Rq is the current curvature Rq ~ (lO _33 eV) 2 . Note that Rq > R due to the contribution from matter. In fact, 
if we can regard f(Ro) as an effective cosmological constant, the effective Einstein equation gives 

Ro = Ro — K 2 T mat ter ■ (10) 

Here T ma ttor is the trace of the matter energy-momentum tensor. We should note that f'(Ro) need not vanish exactly. 
Since we are considering the time scale of one-ten billion years, we only require \ f'(Ro)\ <C (lO _33 eV) . The last 
condition corresponding to the second one in ([5]) is: 

Hm/(i?)=0. (11) 

In the above class of models, the universe starts from the inflation driven by the effective cosmological constant ([5]) 
at the early stage, where curvature is very large. As curvature becomes smaller, the effective cosmological constant 
also becomes smaller. After that radiation/matter dominates. When the density of the radiation and matter becomes 
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smaller and the curvature goes to the value Rq ©, there appears a small effective cosmological constant (J9)). Hence, 
the current cosmic expansion could start. 

Eq.Q indicates that there could appear an anti-gravity regime when 1 + f'(A) = 1 + f'(R) < 0. If we assume 
that the anti-gravity does not appear through the known universe history (even in future), the condition f'(R) > — 1, 
combined with the condition (fTTj) . gives 

f(R) > -R , (12) 

which could contradict with ([9]). Since 70% of the total energy density of the present universe could be dark energy, 
we find i?o ~ 0.7-Ro, that is 

f(R ) ~ -1AR , (13) 

which seems to conflict with Eq. (TT2"|) . The problem could occur even for the model [9]. Note, however, (fT2"|) does 
not always mean that the anti-gravity regime occurs at current universe. The condition (|11[) might be a condition 
in future and therefore the anti-gravity might appear in future when the curvature becomes smaller than the present 
value. More detailed discussion will be given later. 

Since the model corresponding to © has been investigated in [TTj ] . we now propose a model which satisfies 



lim f(R) = aR m , (14) 



R 

with a positive integer m and a constant a. The condition to avoid the anti-gravity f'(R) > — 1 tells a > and 
therefore f(R) should be positive at the early universe. On the other hand, Eq.© or (fT3"|) shows that f(R) is negative 
at the present universe. Therefore f(R) should cross zero in the past. 

At the early universe, if the scalar curvature is large, the /(i?)-term could behave as (fT4|) and would dominate the 
Einstein-Hilbert term (if m is bigger than 1). Now let assume that there exists matter with an equation of state 
parameter w. For spatially flat FRW universe, 

ds 2 = -dt 2 + a(t) 2 ( dxl ) 2 - ( 15 ) 

i=l,2,3 

as shown in Q, the scale factor a(t) behaves as 

a(t) ex t h ° , h = 2 "\ . (16) 
3{w + 1) 

Then the effective equation of the state parameter w c g, which is defined by 

Wcff = -1 + ^- , (17) 

can be less than —1/3 and the accelerating expansion could occur if m is large enough although w Q g > — 1. Then the 
inflation could occur due to the R m behavior of f(R) in (fl"4"|) . The expansion of the universe is quintessence like since 
w e ft > — 1. We need not, however, a real quintessence field as the inflaton. Note, w can vanish, which corresponds 
to dust, that is, cold dark matter or baryons, or w can be equal to 1/3, which corresponds to the radiation. If m is 
chosen to be large enough, ho also becomes large and therefore u> e ff goes to — 1 , which coresponds to the cosmological 
constant. Hence, for such class of models, at the early stage, the universe starts from the inflation driven by the 
i? m -behavior, instead of the effective cosmological constant in (J5J. Since w c ff > —1, the curvature could become 
smaller as time passes and f(R) could cross zero. At that stage, f(R) could be neglected and the radiation/matter 
could dominate. After that, f(R) becomes negative and its absolute value could increase as the curvature becomes 
smaller. When the density of the radiation and matter becomes smaller due to the expansion of the universe and the 
curvature goes to the value Ro ©, there could appear a small effective cosmological constant ©. Hence, the current 
cosmic expansion starts. Note that another possibility to realize the above scenario is to add sub-leading curvature 
terms (dominant over R) in such a way that de Sitter like inflation occurs. Moreover, the reconstruction of such f(R) 
gravity at large curvature may be done so that inflationary stage is instable and curvature-induced exit occurs. As 
the last possibility to achieve the exit from the inflationary era one can add (small) non-local gravity action of the 
sort recently proposed in refs. [t3. [l3l| . 

A simplest example satisfying the above conditions is 

f(R) = aR m - (3R n . (18) 
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Here a and j3 are positive constant and to and n are positive intergers satisfying the condition m > n. Since 

f'(R) = amR™- 1 ~ flnFC 1 " 1 , (19) 

we find 

Re= — , f(Ro) = -fj(l--)[^ L ) <0. (20) 

Since 

f(R ) = -2R ~ R ~ (1(T 33 eV) 2 , (21) 

one gets a ~ Ro~ m an d /3 ~ ^o~™- This shows that /(-R) becomes larger than R, f{R) S> i?, even in the solar system 
where R ~ 10~ 61 eV 2 , which could be inconsistent. 
Another, more realistic proposal is 

a R m + l - RR n 

ftR) " I+W ■ (22) 

Here a, /3, and 7 are positive constants and to, n, and Z are positive integers satisfying the condition m + I > n. (If 
necessary, to achieve the exit from the inflation more curvature terms with powers less than to + I may be added to 
Eq.([22]).) For simplicity, we now choose 



Then since 



one gets Rq satisfying is given by 



to = I = n . (23) 



f (R) = K — * - 1 , 24 

(1 + 7-R") 



l/n 



7/ \ V a 



^ = <l: + , (25) 



and therefore 



/(iZo) ~ ~2R a = 4 I 1 + " ^ I • (26) 



As a working hypothesis, we assume Pj/a 3> 1, which will be justified later. Then we have 

Ro ~ (—) ' , f(Ro) = ~2Ro ~ -- ■ (27) 

One also assumes 

f(Rj) ~ ^ j Ri~R T , (28) 

when i? is given by the scale of the inflation, i? ~ i?/, which is Ri ~ (lO 15 GcV) 2 = (l0 24 eV) 2 . The above conditions 
([27)) and (JH]) could be solved as 

a ~ 2R a R„ 2n , /3 - 4Hg J Rj 2n iZy _x , 7 - 2R a R Q 2n R™- 1 . (29) 
Then we find £7/ a - iR^R^ 2n i?| n " 2 - Kpsfr-i^ which is surely large if „ > 1. 
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The action ([3]) may be presented in scalar-tensor form. By using the scale transformation g M „ — > e a g^ v with 
<7 = — In (1 + f'(Aj), the Einstein frame action follows Q: 



V{a) = e°g (e"') - e 2 */ (g ( e -)) = ^ - . (30) 

Here = i? + /(#) and g (e~ (T ) is given by solving a = - In (1 + /'(A)) = lnF'(A) as A = g (e^). After the 

scale transformation g^ v — > e a g^, there appears a coupling of the scalar field a with the matter. The mass of a is 
defined by 

a _ 1 rf»y(°Q _ 1 / A 4F(A) i 1 

ff ~2 rfd 2 2|F'(A) (F'(A)) 2 F"(A)J ' 1 j 

Unless rn CT could not be large, there appears the large correction to the Newton law. 

We now investigate the above correction to the Newton law in the model (|2"2"]l with (|2T)|) . In the solar system, where 

R~ (l(T 61 eV)' ', or in the air on earth, where R ~ (10 50 eV) , we find 

F(R) = R + f(R) ~ R-2R ~ R , 

F'(R) = 1 + f'(R) ~ 1 + — i?"- 1 ~ 1 + n (4- 

7 \Ri 

n(n — 1) / R N 



F"(R) = f"(R) ~ -^-^ ^— j . (32) 

Then in the solar system, the mass m a of the scalar field a is given by m 2 ~ iQ-i60+i09n e y2 anc j j n a - r on 
earth, to 2 ~ l[p 144 + 98n eV 2 . In both cases, the mass m a is very large if n > 2 and the correction to the Newton law 
is very small. 

Thus, it is proposed a viable modified gravity which unifies curvature-induced R m inflation with effective Lambda- 
CDM cosmology. There is no violation of Newton law in such theory while the known universe expansion history is 
reproduced. 



III. ACDM ERA AND ITS FUTURE EVOLUTION 

Let us discuss further properties of the pro pos ed modified gravity. There may exist another type of instability 
(so-called matter instability) in f(R) gravity \l4 . Il5j . It is known that it is absent in the model of ref.Q. The 
instability might occur when the curvature is rather large, as in the planet, compared with the average curvature at 
the universe R ~ (lO" 33 eV) . By multipling Eq.® with g» v one obtains 

0-HR I F<3) ^) y r^r i F '^ R 2F ^ K " t m) 

°- DR+ WW) WpRW R+ 3F^W)~ ^JORT 6F^{R) T ■ (33) 

Here T = and F^(R) = d n F(R)/dR n . We consider a perturbation from the following solution of the Einstein 

gravity: 

R = R b = --T > . (34) 

Note that T is negative since |p| <C p in the earth and T = —p + 3p ~ —p. Then we assume 

R = R b + R p , (|i? p | < \R b \) . (35) 

Now one can get 

= -dfRp + U(R b )R p + const. , 



u{Rb) = {fIvIrT)- F(i)(R b y) VpRbVPRb ' 3 

F^(R b )F^(R b )R b F^(R b ) , 2F(R b )F^(R b ) F^(R b )R 



I: 



3F( 2 )(J? fe ) 2 3FW(Rb) 3F( 2 ){R b ) 2 3FW(R b ) 2 



(36) 
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Then if U(Rb) is positive, R p becomes exponentially large as a function of t: R p ~ ev^'' and the system becomes 
unstable. In the model (22) with (29). if n > 2 

Therefore there is no such instability in the model under consideration. 

Let us consider what occurs when f'(R) — > — 1. If f'(R) < — 1, the theory enters anti-gravity regime as is seen in 
(3). In the (effective) FRW equations with flat spatial part, 

^T = p, = -^(2H + 3H 2 )+p, (38) 

anti-gravity means negative k 2 s . When k 2 s < 0, there is no solution of the FRW equation (58) . which means that 
the anti-gravity could not occur in the FRW universe with flat spacial part. Now we assume f'(R) = — 1 when 
R = R A > 0. When f'(R) -> -1 but f'(R) > -1, it follows K^g -> 0. Then from (38), i?, i? -> when rc 2 ff -► 0, 
which seems to contradict assumption Ra > 0, since the scalar curvature R vanishes when H = H — 0. This indicates 
that the scalar curvature R could not reach Ra- Then, anyway, the anti-gravity could not be realized for the real 
universe even in the future. A possibility of the transition between normal gravity and anti-gravity might be if p and 
p vanishes when k c s vanishes, then H and/or H might not vanish. As usual matter gives positive contribution to p, 
one needs the negative contributions to p. One contribution might come from the negative cosmological constant and 
another might come from the negative spatial curvature, which gives a contribution to p as — 1/a 2 . There could be 
one more possibility for the transition between normal gravity and anti-gravity, where H vanishes but H is finite, and 
therefore the scalar curvature does not vanish. We should note that p could be positive in the flat spatial geometry 
but p can vanish or even can be negative as for dark energy. Hence, one can speculate that the pre-inflationary era 
may result from the transition from anti-gravity to usual f(R) gravity at the point with zero effective Newton coupling 
and infinite negative cosmological constant. 

It is interesting to investigate if one can distinguish the ACDM epoch from usual GR and the same epoch which 
appears in the present /(-R)-model. The analog of the first FRW equation is: 







■^p- + 3 (H 2 + ij) F'(R) - 18 (4H 2 H + HH^j F"{R) + K 2 p mattcl . . (39) 

For the constant equation of state matter, it is known that /9 ma ttcr = paa^ 2 /^ w+1 \ Proposing ([5]), f(R) can be 
expanded with respect to R — Rq as 

f(R) = -2Ro+ Sf , 6f = f (R-Ro) 2 + o((R -R ) 3 ) . (40) 



Here fo is a positive constant. If we keep only the first term in (|40p by putting Sf — 0, we find the following solution 
in ((391) 

« . , «,) . „(„ s (asw (5^)) ■ (4D 

Here A 2 = poa Q 3 ^ 1+w ^ / Rq, I 2 = 3/i?o- The solution (|4"Tj) is the same as for Einstein gravity with a cosmolgical 
constant and matter. One now treats Sf as a perturbation, which could be justified in near future or near past. By 
putting g(t) = ga(t) + Sg and using (|4T3|) , Eq. (|3"9"]) gives 

= -6H Sg - 2l f f Pm ^!° Sg-hf + 3 (ff 2 + H ) Sf 
3(w + 1) 2 V / 

-m^Ho + HoHoJSf" + 0(Sg 2 ) . (42) 
Here subindex "0" expresses a quantity given when Sf = 0, especially 

Ho = 90 (t) , /Wtero = Po% W+^e-W+m® . (43) 

The solution of (32) is 

sg - e-^r*'**!^ f^ ) ( c ls m + (3H (tr + M^)sm 

-18 (flo(t) 2 tfo(f ) + H (t')H (t'j) Sf"(t')) e^FTj f + q ( d f) . (44) 
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Let the present time t = to- When t ~ to, we may assume Hq, Hq, and Hq are constants. Furtheremore one may put 
Sf oc (R — R ) 2 ^ and 5/' oc i? — i? ~ 0, and (5/ = 2/ . Then at leading order with respect to t — to, Eq. (|H)) has 
the following form: 



6 UH Q H + Ho) fo 

5g go (t - t ) , 5o = hJ(T+w) ' ^ 

which gives the correction to the standard ACDM model from /(i?)-gravity. Here it is assumed Sg = when t = to- 
For the model with and (|2l?]). one finds 

fo = °" a ^- 2 (7flS-D . (46) 
(l+ 7 i? f 

From dHJ), we find jR^ - 2RoR^ n Rj^ 1 > 1. Then @6J| could be approximated by 

an 2 R n ~ 2 

fo ~ . (47) 

7 

Then, the order of go may be estimated as 

go = O (r^^/R^- 1 ) ~ l(r 114 " +81 eV . (48) 

Especially for n = 2, we find 50 ~ 10 _147 eV. Since 10 Gigayears correspond to (lO -33 eV) , Sg will be of the order 
of unity only 10 115 Gigayears later. Then the correction does not seem to be observable in the near future or past. Of 
course, the linear approximation used in (|42|) or (|45|) may be not enough in the model (|22|) with ([29]) and non-linear 
terms account may be necessary. 

In the search for footprints of non-linear modified gravity at the current epoch one should not forget that the 
action of such model may be further modified in the future so that its good properties survive. One possibility is the 
following term in the action: 

Sf(R) = . (49) 

Here p and q are positive integer and rj and C are positive constants satisfying the conditions £ <C Ro +q and 77 -C Ro +1 ■ 
We should note that 6f(R) satisfies the condition pip . Hence, at the present universe 

\Sf(R )\ ^ryi?-««i? , (50) 

and, therefore, the <5/(i?)-term can be neglected. When the curvature is larger than Ro, \ f(R)\ is smaller, which shows 
that the <5/(i?)-term is irrelevant in the past universe. (Note, however, it may be used to improve the observational 
predictions of modified gravity when it is necessary). If the curvature becomes smaller and satisfies the condition 

C i/(p+«) « R « Rq > ( 51 ) 
Sf(R) ~ -7?ir« . (52) 

R q o +1 » V » RoR q » RoC q/(p+q) (53) 

Sf(R) could dominate at the future universe. Using the arguments of ref.[J], if matter is included, we find the future 
universe may enter to phantom era: 

«(«)-(*.-*)*•, (54) 

In the above a(i), there seems to appear the Big Rip singularity at t — t s . Near the singularity, however, the curvature 
becomes large and (5/(i?)-term does not dominate. Therefore, the Big Rip singlarity does not occur and the phantom 
era is transient one. After that the universe enters a quintessence or ACDM epoch. 

Our consideration shows that even if the current universe (as predicted by modified gravity under consideration) is 
qualitatively/quantatively the same as the standardACDM era this may not be true in the near future. For instance, 
the (transient) phantom epoch may emerge without the need to introduce phantom matter. 



5R behaves as 



If rjR q ^> Rq, which requires, 
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IV. DISCUSSION 



In summary, we proposed a modified f(R) gravity which predicts natural unification of early-time inflation with 
late-time acceleration. This theory which is closely related with the models [ll[ passes the local tests (Newton law, 
stability of Earth-like gravitational solution, heavy mass for additional scalar degree of freedom, etc). The speculative 
possibility of past or future anti-gravity regime is briefly mentioned. The evolution of the f(R) ACDM epoch is 
discussed. It is shown that it may get out from the cosmological constant boundary by future reconstruction of the 
gravitational action. As a result, the future universe may enter a quintessence-like or transient phantom era or it 
may continue to be asymptotically de Sitter universe forever. At the next step, it is necessary to investigate the 
cosmological perturbations in the highly non-linear gravity under discussion. However, this is a quite non-trivial task 
as the results should be presented in a gauge-independent formulation (the simple approximation which is analogous 
to the one made in GR with dark fluid does not lead to realistic predictions due to its gauge dependence). As 
more precise observational data for cosmological parameters are expected very soon, the further study of various 
cosmological predictions of our model (to distinguish it from GR) are requested. This will be done elsewhere. 
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APPENDIX A: CONSISTENCY OF STELLAR SOLUTION IN /(E)-GRAVITY 

Let us investigate if /(i?)-gravity admits a consistent stellar solution, where vacuum solution matches onto the 
stellar-interior solution. 

For this purpose, one may rewrite Eq.([2]) in the following form 

\g^R - R„ v - ~<?^A + yT (m) ^ = - X - 9vuV (f(R) + A) + R^f'iR) + g^af(R) - V^„f'(R) . (Al) 

Here A is the value of f(R) in the present universe, which corresponds to the effective cosmological constant: A = 

/(i?o), then we have Rq ~ A ~ (l0 _33 eV) 2 . We now like to treat the r.h.s. part as a perturbation. The last two 
derivative terms could be dangerous when we consider the stellar configuration since there could be a jump in the 
value of R on the surface of the star. One may regard the order of the derivative could be the order of the inverse of 
the Compton length of a typical scale of the system. Since the most dangerous case corresponds to particle, one may 
estimate the order of the derivative could be the Compton length of proton: 9 M ~ m p ~ 1 GeV ~ 10 9 cV. Here m p is 
the mass of proton. We also assume the scalar curvature has the order of R e ~ 10 _47 eV 2 , which corresponds to the 
curvature inside the earth. 

First we consider the 1/R model: 

f(R) = ~ ■ (A2) 
The order of the dimensional parameter fi is 10~ 33 eV. One may estimate 

nf(R) ~ V„V v f{R) ~ ^ ~ 10- 20 eV 2 , (A3) 

which is much larger than A or R e and therefore the perturbative expansion breaks. 
In case of the HS model ([3]), we find 

Df'(R) ~ V M V„/'(ii) ~ ^ (A;) ~ io-3-i7» eV 2 . (A4) 

In (g}, R/m 2 — 41. Eq. (|A4p shows that if n > 2, Of(R) or V M V„/(.R) could be much smaller than R and therefore 
the perturbative expansion is consistent. For the model (|22[) the qualitative structure is similar to that of the HS 
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model. This shows that for the class of models under discussion, the stellar solution is qualitatively similar to the one 
in Einstein gravity. 
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